Controlling a mechanical oscillator with a tunable coherent feedback network 
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We demonstrate a fully cryogenic microwave feedback network composed of distinct superconduct- 
ing devices interconnected by waveguides and designed to control a mechanical oscillator coupled 
to one of the devices. The network is partitioned into an electromechanical device to be controlled 
and a dynamically tunable controller that coherently receives, processes and feeds back continuous 
microwave signals that modify the dynamics and readout of the mechanical state. While previous 
electromechanical systems represent some compromise between efficient control and efficient readout 
of the mechanical state, as set by the electromagnetic decay rate, this flexible controller yields a 
closed-loop network that can be dynamically and continuously tuned between both extremes much 
faster than the mechanical response time. We demonstrate that the microwave decay rate may be 
modulated by at least a factor of 10 at a rate greater than 10 4 times the mechanical response rate. 



INTRODUCTION 

Feedback methodologies underlie or are a key to under- 
standing robust systems in all fields of science and tech- 
nology. Analyzing complex natural systems (e.g. bio- 
logical processes, financial markets, Earth's climate) and 
engineered systems (e.g. electrical circuits, autonomous 
vehicles, utility grids) in terms of real-time feedback il- 
luminates their dynamics. Moreover, the proper applica- 
tion of feedback is an essential ingredient in the successes 
of James Watt, the Wright brothers, automated indus- 
trial production, scanning tunneling microscopes, and a 
decentralized internet, to name a few pQ. Researchers are 
rapidly improving their ability to engineer fully coherent 
electromagnetic (EM) devices that operate at the quan- 
tum level, but given the extreme demands of quantum- 
scale engineering one must anticipate that feedback con- 
trol will be essential in realizing useful complex systems 
[2H1] . Unfortunately, despite some progress [4UT6] , a sys- 
tematic craft of control design using coherent systems is 
lacking. Coherent EM control design is more restrictive 
than traditional electrical circuit design since the tech- 
nologies used must be describable by unitary processes 
[SHE]- As such, fully coherent EM control systems de- 
signed to tune even classical dynamical processes [TTJ [13] 
can help advance researchers towards the goal of a sys- 
tematic craft of quantum control. 

Here, we demonstrate a feedback network of coherent 
EM devices: a coherent superconducting electromechan- 
ical device [17-24 to be controlled (the "plant") and a 
coherent superconducting microwave controller [13] [25] 
with which it exchanges continuous coherent signals and 
is thus regulated. This coherent control approach is 
roughly analogous to using Watt's centrifugal governor 
to regulate a steam engine, in that the controller and 
the plant are built using similar technologies and the 
dual processes of sensing and actuation are effectively 
inseparable [TTHT3] . Such an approach is distinct from 



measurement-based approaches, in which coherent sig- 
nals produced by the plant are measured, the measure- 
ments are interpreted by an incoherent controller (often 
a computer), and the controller synthesizes new signals 
that act back on the plant [26-32 . In comparison, coher- 
ent control systems are more physically compact. They 
also make better use of given resources — e.g., a coherent 
controller has access to both quadratures of a coherent 
signal, without any added noise [3 [9]. 

The field of electromechanics (also known as optome- 
chanics) concerns mechanical oscillators coupled to itin- 
erant EM fields [33]. While our network is too complex 
for traditional electromechanical models [33 -35 , it re- 
quires only three components, all of which are accessible 
to superconducting circuit labs, and is efficiently and in- 
tuitively modeled as a linear feedback system [TJ [5j [8] . 
This network was inspired by a new, coherent feed- 
back interpretation of a common electromechanical tech- 
nique [13 HU [22] [33] ESJgJj], but it provides a previ- 
ously unavailable degree of dynamic flexibility in how an 
electromechanical oscillator may be controlled and mea- 
sured. For instance, while previous electromechanical 
systems make a compromise between better mechanical 
control and better mechanical measurement capabilities, 
we demonstrate that this network may be dynamically 
modulated between both extremes. This is equivalent to 
alternately isolating the oscillator from and coupling it 
to fields outside the network, without modifying its cou- 
pling to intra-network fields. In future work, the wave- 
form of signals entering or exiting the network could be 
modulated, facilitating electromechanical state transfer 

mum. 

Sideband cooling is a widely used electromechanical 
technique [T7] EE] [22] [32 [36HI0]. The canonical side- 
band cooling system consists of an EM resonator cou- 
pled to a mechanical oscillator and driven at a frequency 
below the EM resonance. Mechanical cooling via this 
method is usually explained in analogy to early ion trap- 
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ping experiments [33H351 145] . A small portion of the 
low-entropy EM drive is inelastically scattered by the me- 
chanical oscillator in the resonator, with each scattered 
photon having gained or lost a phonon's worth of energy. 
If the center frequency of the resonator is higher than the 
drive frequency, then more frequency up-converted pho- 
tons are scattered than down-converted photons. Con- 
sequently, the oscillator gives up energy on average to 
the EM field per scattering event, reducing its motion. 
This phenomenon is usually visible through "blue" (up- 
converted) and "red" (down-converted) sidebands to the 
"coupling tone" (drive carrier) emitted by the resonator. 
A cooling oscillator is indicated by more blue than red 
sideband power, less total sideband power, and sideband 
power spectral densities broader than the intrinsic oscil- 
lator linewidth. 

There are also coherent feedback interpretations of the 
same system |E1 H6j [47] , but they are seldom emphasized. 
For instance, one may view the input coupler to the res- 
onator as a controller that receives a continuous stream 
of information about the mechanical motion and synthe- 
sizes a reactive feedback signal counteracting this mo- 
tion. If the resonator is a single-sided Fabry-Perot cav- 
ity and the oscillator is the spring-mounted mass of the 
perfectly reflecting mirror [28j [33j [33 [38], for example, 
the partially-transmissive input mirror acts as a coherent 
controller. As the mirror mass moves forward or back- 
ward, optical signals reflecting off it are Doppler shifted 
and arrive at the input mirror with a positive or negative 
phase shift, respectively. The input mirror is designed 
such that, through interference with an external drive, 
more or less light is coupled into the cavity when the 
phase shift is positive or negative, respectively. Higher 
intracavity power counteracts forward motion of the mass 
through radiation pressure, while lower intracavity power 
encourages a retreating mass to be pushed forward by 
its spring. Thus the mass feels a light-induced drag, re- 
ducing its motion. The analogous, measurement-based 
control approach would have no input mirror, but the 
phase of light reflecting off the mass would be measured 
and used to apply a counteracting force to the oscillator 
[7J [9j |28]. To our knowledge, it is novel to identify the 
input coupler as the controller, but with this perspective 
it is clear that employing a single, partially transmissive 
mirror is but one control design. Instead, the controller 
could be an entire coherent subsystem [7J [9] , offering dif- 
ferent types of control over the oscillator. 

The choice of an input coupler in an electromechani- 
cal system typically represents a compromise. A weaker 
coupler (e.g. more reflective input mirror) builds up a 
higher circulating power inside the resonator per Watt 
applied externally. It thus produces an electromechan- 
ical system capable of high feedback "gain" with rela- 
tively low applied power [2TJ [40] (is efficient at "con- 
trol"). Resonators with weak couplers also effectively 
"integrate" Doppler-shifted signals for longer, and can 



consequently respond to finer mechanical motion. How- 
ever, weak couplers frustrate high-bandwidth detection 
of mechanical motion by low-pass filtering mechanical 
information on EM signals leaving the system and en- 
couraging loss of the same through parasitic channels (is 
less suited for "measurement"). Conversely, relatively 
strong couplers produce systems that favor mechanical 
measurement over control [28j [48] . A "weak" ( "strong" ) 
coupler is one that produces a resonator linewidth n t nar- 
rower (broader) than the center frequency of the oscilla- 
tor yielding what is known as a "resolved-sideband" 
("unresolved-sideband") system. In both choices, the 
coupler-controller is both passive and (in all electrome- 
chanical systems to our knowledge) fixed at the time of 
construction. Here, we take an alternate approach and 
control the oscillator with a dynamically tunable coher- 
ent subsystem, enabling a wider range of control and 
measurement capabilities. 

Generically, many superconducting microwave realiza- 
tions of electromechanical systems consist of a high qual- 
ity microwave LC resonator, in which the capacitance is 
modulated by a spring-loaded mass [T71 [I9j [2TH241 148] . 
Motion in this mechanical oscillator adjusts the capaci- 
tance of the resonator, and thus acts on microwave pho- 
tons in the resonator. And the voltage potential across 
the capacitor creates a Coulomb force that acts on the 
oscillator, and thus the number of microwave photons in 
the resonator acts on the position of the oscillator. Via an 
inductive transformer, the intrinsically low-loss resonator 
couples weakly to a 50Q transmission line, which serves 
as the input-output (I/O) port through which the elec- 
tromechanical device is probed, driven and controlled. 
This input transformer, through which microwave fields 
internal and external to the resonator interfere, tradition- 
ally serves the same role as the input mirror in the Fabry- 
Perot system described above. This microwave system 
is accurately described at a quantum level with a model 
that is formally equivalent to the models used to describe 
electromechanical systems of all scales and realizations 
[T71 [34| [35j [39]. Although the specific electromechani- 
cal circuit (EMC) used here is somewhat incidental, it is 
a new design of this type [24 , built from a 4.672 GHz 
center frequency and over-coupled n t /2i: = 2.8 MHz res- 
onator, employing an oscillator with an effective mass 
of ~10 ng, center frequency of Q = 2irx 713.6 kHz and 
intrinsic linewidth of To = 27rx0.81 Hz, and in which 
each photon exerts a force on the oscillator (in coherent 
amplitude units) of go — h x 2.3 Hz. 

Many superconducting microwave electromechanical 
experiments conducted by us [T7J [2TJ [4lJ [48] and others 
HH [19j [22j [39] are partitioned into an "upstream" elec- 
tromechanical circuit (EMC), containing an oscillator to 
be studied, and a "downstream" cryogenic low- or near 
quantum limited-noise temperature amplifier, for high- 
fidelity readout of microwave probes. Our group typ- 
ically uses near quantum-limited Josephson parametric 
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amplifiers (JPAs) for readout [25]. JPAs are composed 
of a 20dB directional coupler followed by a single-port 
"tunable Kerr circuit" (TKC) [13] [25] . a nonlinear mi- 
crowave resonator whose center frequency is drive power- 
dependent and tunable with an applied magnetic flux. 

Housed in a dilution refrigerator, our network is the 
interconnection of a single-port EMC, a JPA, and a mi- 
crowave tee (Fig, [TJl). The EMC and JPA are connected 
to one tee port each, while the remaining tee port serves 
as the overall I/O network port. The network output is 
amplified by a cryogenic HEMT for analysis. While the 
EMC and TKC were originally intended for different ex- 
periments and were mounted in separate sample boxes, 
and interconnections were made using ~cm-length copla- 
nar waveguides and low- loss SMA cables, an analogous 
network could be fabricated on a single Si substrate. In 
this report, we do not pump our JPA; effectively, it is 
operated as a gain-1 amplifier, or more precisely, as a 
linear resonator whose center frequency can be tuned by 
applying either magnetic flux or a moderate amount of 
microwave drive power. We do, however, typically em- 
ploy a cancellation tone that prevents the carrier of any 
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FIG. 1: a) An experimental schematic of the coherent feed- 
back network. A microwave tee connects an electromechanical 
circuit (EMC) to a Josephson parametric amplifier (JPA) — 
here broken down into its directional coupler (DC) and tun- 
able Kerr circuit (TKC) subcomponents — and provides an 
overall input-output port. The In port is typically driven by 
a strong electromechanical coupling tone, but this tone is typ- 
ically cancelled at the DC, preventing its strong carrier from 
driving the TKC nonlinear. Microwave signals are detected 
at the Out port, b) Conceptual schematic of coherent inter- 
connections. Ideally, the JPA acts as an over-coupled and 
wide-band linear resonator with a tunable center frequency. 



strong, electromechanical coupling tone driving the net- 
work from reaching the TKC and driving it nonlinear. 

Conceptually, the network's operation is depicted in 
Fig. [TJd. Microwave signals leaving the EMC and car- 
rying information about the motion of the oscillator are 
split three ways at the tee: most of the amplitude is split 
evenly between the network output and the JPA input, 
and a small amount is reflected. Modeling the JPA as 
a linear resonator with a tunable center frequency, the 
JPA reflects its signal portion with a tunable phase shift. 
Similarly split by the tee, a portion of the JPA-reflected 
signal interferes with the EMC-to-output signal, enhanc- 
ing or diminishing the rate of mechanical information 
leaving the network. The portion of the JPA-reflected 
signal fed back to the EMC can exert a force on the me- 
chanical oscillator, amplifying or counteracting its mo- 
tion. Both effects are controlled by a single, continuous 
parameter — the JPA center frequency — and both ef- 
fects improve as the JPA bandwidth exceeds the band- 
width of the EMC and the EMC bandwidth exceeds the 
oscillator's center frequency. This network is thus anal- 
ogous to the Fabry-Perot optomechanical system, but 
with a user-controlled knob that simultaneously modi- 
fies the reflectivity and position of the input mirror [42] . 
This network is also naturally interpreted as a feedback 
control system: rather than make a room-temperature 
measurement of the JPA output, this signal is fed back 
directly to the EMC without leaving the cryostat. Ide- 
ally, the JPA's transformation of the microwave signal is 
unitary (i.e. adds no noise), and this feedback network is 
fully coherent and thus fundamentally different from any 
measurement-based control system [5H2]. 

MODEL 

Each subcomponent is separably described by an I/O 
dynamical model, a familiar concept in quantum optics 
[5j |49] and superconducting microwave systems [13] [50] . 
Moreover, despite fundamental nonlinear it ies, the de- 
vice dynamics essential to the network's operation are 
accurately described by common linear approximations 
[34] 135] , as discussed more below. In an I/O formula- 
tion, the tee is a 6-input and 6-output (6-1/0) itiner- 
ant field scattering device (three physical ports, each of 
which is both a field input and output, each field de- 
scribed by two quadrature degrees of freedom) with no 
internal degrees of freedom. The JPA is modeled as a 
2-1/0, over-coupled linear resonator device (whose cen- 
ter frequency is a tunable parameter) with two internal 
degrees of freedom, representing the two quadratures of 
the mode inside the TKC [50]. And the EMC is 
a 4-1/0 device, with two I/O ports connected to a mi- 
crowave field (effectively at zero temperature) and two 
I/O ports connected to a thermal mechanical bath. The 
EMC also contains four internal degrees of freedom, rep- 
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resenting the quadratures of coupled microwave and me- 
chanical modes [9] [46] [47] . Fundamentally, when vacuum 
fluctuations in the microwave fields are significant, each 
instance of microwave dissipation inside each device or 
interconnection should be modeled using an additional 
I/O port pair [5] (and a beamsplitter component inter- 
rupting lossy interconnections), but such accuracy is not 
required below. It is well known that a network of linear, 
coherent I/O devices coupled to itinerant Gaussian fields 
may be modeled as a classical control system of inter- 
connected, linear state space models [6H9]. The residual 
"quantumness" of these systems is captured by nonclas- 
sical "noise" driving the network inputs, but again, such 
quantum-level accuracy is not required below as the dy- 
namics considered are classical. A linear network model 
representing the schematic depicted in Fig. consists 
of ten coupled linear equations of motion, with six in- 
ternal degrees of freedom, and two microwave and two 
mechanical bath I/O port pairs. Despite this complex- 
ity, standard control systems software toolkits make the 
construction of network dynamical models intuitive and 
efficient linear system theories may be applied to their 
analysis. 



While the electromechanical interaction is fundamen- 
tally nonlinear, the dynamics we consider may be mod- 
eled by linearizing the coupling of the microwave and me- 
chanical modes about a large microwave coherent state 
induced in the resonator by the coupling tone [34] [35] , 
an approximation that is nearly ubiquitous in electrome- 
chanics. In a frame in which the EM degrees of freedom 
are rotating with the carrier frequency of the coupling 
tone, the Hamiltonian that describes the internal dynam- 
ics of the EMC in the linearized approximation is (h = 1) 

H M = £la\ai + A<4a 2 + g(ai + a[)(a 2 + a\) (1) 

where a± (<22) is the mechanical (microwave) mode an- 
nihilation operator (a[ 2 the creation operator), Q is the 
center frequency of the oscillator, A = uj r — uj c is the fre- 
quency detuning between the resonator center frequency 
(uo r ) and coupling carrier (cj c ), and g = g§^Jn~ c , where n c 
is the number of intra-resonator photons induced by the 
coupling tone. 



In the usual I/O formulation for an inductively coupled 
EMC, the Heisenberg equations of motion may be repre- 
sented in a rather condensed form by virtue of the fact 
that they describe a dynamical system that is both linear 
and fundamentally unitary [HJ [9] . These restrictions on 
form ultimately stem from the need to preserve canonical 
commutation relations and are the primary distinction 
between coherent and more general linear I/O systems. 
For instance, we may write the I/O equations of motion 



for the EMC as 
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dt 
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Am cl + Bm bi n 
Cm o> + Dm b in 
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where a 



<2i, (22, cl\, a\ 



(T indicating transpose), 

T 

, and similarly for b out . 

Here, b T h,in and b E ,in (p T h,out and b EjOU t) are standard 
I/O theory objects that are functions of time t and may 
be roughly considered annihilation operators on the in- 
finitesimal segment of free field incident on (leaving from) 
the device at time t [8 . Specifically, b Th ^ in / out and 
bE,in/out are associated with the fields that drive the sys- 
tem through the thermal mechanical port and the EMC's 
microwave port, respectively. The matrices Am, Bm, 
Cm, and Dm define the dynamics and despite their size 
are themselves represented by a relatively small number 
of variables. For instance, for general linear coherent sys- 
tems in which the matrices C and D are block diagonal 
(i.e. no b out is an explicit function of any b\ n or a)) these 
matrices may be broken down as [8] 
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Am = —^C M Cm + jn M , B M = —C M D M , 



C 



M 



diag C M ,C M , D M = diag D M ,D M (3) 



where * is the element- wise complex conjugation and j 
is the imaginary number with the electrical engineering 
sign convention [51]. The dynamics are thus completely 
defined by three matrices 
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where r is the intrinsic mechanical decay rate, n c is 
the EMC coupling (internal) decay rate. These expres- 
sions exemplify the high symmetry in all linear coherent 
I/O models. They may be reduced even further [8 , but 
it doesn't suit our purposes to do so here. 

Eqs.[2]are of a form equivalent to the state space repre- 
sentation of a linear system common in classical control 
theory. In a control context they would be succinctly ref- 
erenced using an object M that is defined by Am, Bm, 
Cm and Dm and is often notated as |TJ [8] 



(5) 





Bm 


_ Cm 


D m 



(each matrix and its location in the above array implies 
Eqs. |2j. We can similarly generate state space model 
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representations for the JPA and tee devices; call them JJ 
and T, respectively. 

We now come to the problem of constructing the net- 
work dynamical model. It has been previously pointed 
out that the formal equivalence between these quantum 
I/O models and the state space models that appear in 
classical linear control systems theory means that anal- 
ogous procedures for constructing dynamical models in 
classical linear control systems also apply to linear quan- 
tum I/O ones [8]. Through these methods, we may derive 
a new state space model N that describes network as a 
whole. In the remainder of this section, it will be ambigu- 
ous whether we are constructing a quantum or coherent 
classical model. It is only when we neglect the effects of 
quantum fluctuation in the free fields when we compare 
our model to experiment that the approach becomes a 
classical approximation [TTJ [13] . 

The network model we need is depicted schematically 
in Fig. [2j as inspired by Fig. [TJd, where the EMC and 
JPA each exchange input and output signals with two 
ports of the tee. The u ip" and "0" blocks represent the 
phase shifts accumulated just from transmission line de- 
lays around the tee-and-EMC and tee-and-JPA network 
loops, respectively (a more physical model would not 
lump these phase delays only on the input side of the 
tee, but half on each input and output. Up to unimpor- 
tant phase shifts to the network internal variables, the 
two modeling approaches yield identical predictions). We 
thus define our network by defining the state space model 
objects T, M, and JJ, and labeling the fields associated 
with each component's inputs and outputs to represent 
the schematic Fig.[2) For instance the label LL E Q " (the la- 
bel "E*") in Fig.gfis applied to both the M output field 
operator bE,out out ) an d the appropriate field driving 
one of the T inputs. The label "/V is applied only to 
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FIG. 2: Control system schematic of the network depicted in 
Fig. [T]d. The tee, EMC and JPA components are represented 
as multi-input multi-output (MIMO) linear state space mod- 
els that continually exchange signals. The interconnections 
are labeled here such that, e.g. E^ stands for the two free 
field labels E Q and E*, which are explained in the text. Trans- 
mission line delays may be modeled here as static phase shifts 
ip and 0. Standard software toolkits are available that reduce 
the network on the left to a single MIMO state space model N. 
This approach is compatible with both quantum and classical 
investigations. 



the (heretofore unmentioned) b^^ n field that drives the 
appropriate T input representing one of the network's 
microwave inputs, etc. While algorithmic, the construc- 
tion of N is extremely tedious. However, in the limit that 
transmission line delays may be accounted for by static, 
dispersionless phase shifts between the subcomponents, 
software toolboxes are available in, for example, Matlab 
and Mathematica that completely automate this model 
construction. This basic fact is under-recognized by the 
Physics community and deserves explicit emphasis: after 
defining the three state space models and labeling the 
component inputs and outputs according to Fig. |2j the 
entire state space model N may be obtained in Matlab 
using the single command N = connect (T,M, J) . (other 
software that can deal with nonlinear quantum I/O net- 
works has also been developed [T3l [52] ) Using such meth- 
ods, these control systems toolkits determine an effective 
network state space model defined by matrices An (di- 
mensions 6 x 6), B N (6 x 4), C N (4 x 6), and D N (4 x 4), 



N= ^& . (6) 



' A N 




_ Cn 





which are too unwieldy to write out in general here. 

It is useful, though, to consider the dynamical model 
for just For the EMC in isolation, and taking g = 0, 
the equation of motion for this mode is 



d_ 
dt 



&2 = (jA - y ) a 2 - j^b E 



(7) 



where K t = ki + k c . In the network N, taking g = and 
in the limit that the JPA has no internal loss and couples 
much more strongly to transmission lines than the EMC 
does, we find that 



d_ 
dt 



a 2 = (jA',, - y) oa + e-* A +i*> (8) 



where n n = k,i + n' c and 
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AL = A- Kc 



(|) sin(V) - sin(6»' + V) 
4cos2(f) s in 2 (f) W (*¥) 



A = arctan ( cot | — 1 — tan I — 



(9) 



where 0' is the total phase shift acquired by a narrow- 
band signal that makes a round trip from the tee to JPA 
and back (i.e. is determined by the length of that inter- 
connection and the center frequency of the JPA). Thus, 
in this approximation, the EMC resonator response still 
acts like a single- mode resonator, except one that is now 
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driven by the inputs "/4 and has a new detuning from 
the carrier tone (A^,), a new decay rate (« n ), and a new 
phase shift between inputs and outputs (controlled by A) 
that are all controlled by 0' ', which is in turn controlled 
by the JPA center frequency. In general, and especially 
when intra- JPA loss is not negligible, these same quali- 
tative effects still hold, although they lose their quanti- 
tative accuracy, in which case the full network model N 
is needed. 

Below, we compare the network model N's steady state 
output predictions to data. This is most naturally done 
through a Laplace transform of the field and internal vari- 
ables such that, e.g. 

/'OC 

*WnM = / e- st b^ in dt (10) 
Jo 

where b a ^ n (as above) is a function of time. In the 
Laplace domain, it is appropriate to write the network 
equations of motion in terms of susceptibilities and trans- 
fer functions such that [8 

a[s] = (si- A^^Bnbinls] = x[*]&m[s] 
b ou t[s] = (C N (sI - An^Bn + D N ) b in [s] = E[s]b in [s] 

(11) 

where a = [a\ ct2 as a\ a\ a\) T where as is the annihi- 
lation operator on the internal JPA mode, b in / out = 

[b T h,in/out b^in/out ^ThM/out ^in/outV' and 1 is an iden " 

tity matrix of appropriate rank. This representation is 
only valid when the system is Hurwitz stable, i.e. the real 
parts of the eigenvalues oi An are negative. While man- 
ually calculating the multi-input multi-output (MIMO) 
transfer function S[s] would be prohibitively tedious, 
again standard software toolboxes completely automate 
the procedure; in Matlab, given a network state space 
model N, the MIMO transfer function representation is 
obtained using the single command Xi = tf (N). 

The transfer function 5[s] is applied several times 
below to make predictions about the signals emitted 
by the network. For instance, the network's phase re- 
sponse to different frequency microwave probes is sim- 
ulated in Fig. [4J3 by calculating the phase angle of 
the {/i Q ,/ii} matrix element of S[s], angle [E{ Mo/x .}[s]] , 
as s runs along the imaginary axis. In general, the 
microwave-in-to-microwave-out response will be referred 
to as the network's Sn response. Similarly, the absolute 
value squared of the {/i 07 Thi} matrix element of 
p{/i Thi} I ? represents the gain with which ^-energy 
excitations in the thermal bath induce ^-energy excita- 
tions in the microwave network output (i.e. induce blue 
sideband power in the network output), through effec- 
tive yOUt bTh, in + h.c.^j -type interactions. Conversely, 
|^{/i*Thi} [jQ] | 2 represents the gain of thermal bath- 



induced ^-energy de-excitations in the microwave net- 
work output (i.e. induction of red sideband power in the 
network output), through effective (b ayOUt bTh,in + h.c.)- 
type interactions. The thermal-bat h-in-to- microwave- 
out response will sometimes be referred to as the net- 
work's Si 2 response. Similarly, 

arg max |S {/ioThi} [j(tt + SQ)] | 2 (12) 
sn 

may be used to predict microwave-induced shifts in the 
mechanical center frequency [TTJ [39] . Finally, an expres- 
sion related to xb'^L ^ ne it 1 0,0*2} matrix element of 

CnUM-An)-^ (13) 

may be used to predict the rate with which oscillator 
excitations decay out the network blue sideband, poten- 
tially enhancing the total oscillator decay rate, and sim- 
ilar expressions may be used to predict red sideband de- 
cay and losses dissipated internally to the network. Thus, 
the functions 5[s] and x[ s ] are extremely convenient for 
modeling steady state network dynamics. And while the 
matrix of expressions represented by S[s] is difficult for 
a human to parse, for instance, the object S[s] is easily 
manipulated algebraically and computationally. 

RESULTS 

With the JPA far-detuned, the network's Sn phase re- 
sponse (i.e. the phase with which microwave tones are 
reflected by the network — see previous paragraph.) near 
the center frequency of the EMC indicates an overcou- 
pled, ft n /27r = 3 MHz linewidth network resonance at 
4.672 GHz, Fig. [3^i. Driving with a strong coupling tone 
on this network resonance (and protecting the TKC with 
a cancellation tone), the thermal motion of the oscillator 
is visible via red and blue inelastically-scattered network 
output signals, i.e. sidebands. The power spectral den- 
sity of these sidebands (Fig. [3J3) indicate a mechanical 
center frequency of Q = 2tt x 713.6 kHz and intrinsic 
damping rate of To = 2tt x 0.81 Hz. By applying different 
amounts of external magnetic flux to the TKC, the appar- 
ent network resonance frequency and linewidth vary to- 
gether, each by values of order the EMC's coupling decay 
rate (Eqs. [9|. As the network's response to probe tones 
varies, the characteristic mechanical motion inferred from 
the sidebands also varies. In particular, as the network's 
center frequency moves higher and its linewidth narrows, 
the apparent mechanical motion is damped more heav- 
ily (sideband linewidth becomes T > r ), is cooled, and 
its center frequency moves (SQ ^ 0). These effects are 
completely analogous to sideband cooling in traditional 
electromechanical systems [33 . In our system, the origin 
of these effects may be interpreted as sideband cooling, or 
as the JPA acting as a coherent controller, exerting dif- 
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ferent control functions for different amounts of applied 
flux. It is important to note that n c does not change 
as the TKC is tuned. While the network's electronic re- 
sponse to probes and mechanical sidebands varies with 
the TKC state, the coupling tone is canceled before the 
TKC and thus cannot be effected by it. 

The network response may be finely tuned. Probing 
the network over a 2 GHz range, the JPA center fre- 
quency varies with applied flux, Fig. [4^. The other visible 
resonances, the narrow-band EMC and very broad-band 
resonances arising from standing waves resonances be- 
tween the subcomponents are only effected as the JPA 
tunes through them. Focusing on the 10 MHz band 
around the point the JPA tunes through the EMC (« 56 
/iA flux coil bias; Fig. [4)3), we see that as the flux in- 
creases, the network resonance first broadens slightly, 
then abruptly narrows and shifts higher, then becomes 
under-coupled, and finally re-broadens and returns to 
its original center. Device and network parameters may 
be precisely estimated given these 2 GHz and 10MHz 
probe ranges. In particular, we infer an EMC linewidth 
of n c = 2tt x 2.7 MHz and internal loss rate ni = 2tt x .1 
MHz; a JPA linewidth of 7 = 2tt x 50 MHz and in- 
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FIG. 3: a) Microwave probe measurements near the elec- 
tromechanical network resonance for three flux bias settings. 
The JPA is far-detuned at 50 /iA coil bias, but near-detuned 
for the others. The brown arrow represents an electromechan- 
ical coupling carrier tone on-resonant with a 50 /iA-biased 
network. The blue and red arrows represent ±714 kHz elec- 
tromechanical sidebands, b) Employing n c = 1.6 x 10 6 , the 
power spectral density of the blue sideband is shown for the 
three flux settings. 



ternal loss rate of 7/ = 2tt x 5 MHz; and tee-to-EMC 
round trip phase shift of ip = 3.8 rad and tee-to- JPA 
round trip phase phase shift of = 0.22 rad. The tee 
is assumed to be ideal and properly terminated, reflect- 
ing l/9 th of the incident power with a 180° phase shift 
and splitting the remaining power between the other two 
ports. Although difficult to distinguish losses within the 
three over-coupled devices from interconnection losses, a 
loss model dominated by EMC and JPA internal losses 
is more consistent with observations than a model dom- 
inated by interconnection losses. Using these device pa- 
rameters and a flux bias-to- JPA center frequency func- 
tion, a linear control systems model built from intercon- 
nected device state space models using the Matlab Con- 
trol Systems toolbox accurately reproduces the network's 
electronic response as the flux varies, Fig. [4^. 

We now further consider the effects of this network 
on the oscillator dynamics. Driving the network at a 
frequency fixed near resonance when the JPA is far- 
detuned, A' e = 2tt x 13 kHz, and with a power such that 
n c = 1.6 x 10 6 is expected, we measure the linewidth and 
center frequency of the steady state mechanical sidebands 
emitted by the network as the flux is tuned, Figs. [5^i&b. 
Starting from a linewidth of T « To with the TKC far- 
detuned, T reaches a maximum of 2tt x 56 Hz, with a 
simultaneous frequency shift of SQ — — 2ir x 24 Hz, when 



a ) Network angle[Sn] as JPA is Tuned 




Probe Frequency [MHz] Probe Frequency [MHz] 



FIG. 4: a) Network's 4-6 GHz Su phase response to mi- 
crowave probes as flux bias is tuned. JPA resonance shifts 
with flux and couples to EMC (4.7 GHz) and standing wave 
(4.2 & 5.5 GHz) resonances as it passes through them, b) 
Detail of response around the EMC- JPA intersection. Points 
marked with pentagon, circle and triangle indicate network 
electromechanical center frequencies for the 50, 56.2 & 56.5 
/iA bias points depicted in Fig. [3] c) Simulation of response 
detail using a linear control systems model. 
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the flux bias is 56.4 fiA. This corresponds to a JPA center 
frequency 3.4 MHz higher than the coupling tone, and an 
undercoupled network response with a ~ 300 kHz inter- 
nal loss rate. Because the A' e ~ coupling tone has no 
significant effect on the mechanical state when the JPA is 
far-detuned, it is natural to interpret the JPA as "control- 
ling" the oscillator as it is tuned through the EMC. Using 
a control-systems model (here employing standard-issue 
Mathematica toolkits) and no fit parameters, predictions 
for the expected mechanical linewidth and frequency shift 
as a function of flux bias are also included in Figs.[5^&b. 

A natural metric for quantifying the efficiency with 
which the network controls the mechanics is the prob- 
ability that a mechanically scattered photon dissipated 
by the network is a frequency up-converted photon mi- 
nus the probability it is a down- converted photon. This 
"cooling efficiency" (CE) is independent of n c , but is 
coupling-frequency dependent and is positive when the 
network cools the oscillator (and is negative when it am- 
plifies thermal motion). The maximum CE obtainable 
over all coupling frequencies is Q/^Q 2 + ft^/4, where 



70 
60 
50 

S 30 

f-H 

20 
10 




Mechanical Linewidth 
vs. JPA Tunings 



Mechanical Center Frequency 
vs. JPA Tunings 





56 58 
Flux Coil Bias [\i A] 



56 58 60 

Flux Coil Bias [\i A] 













0.5 










® 








-0.5 



"°'l668 4670 4672 4674 
Coupling Tone 
Frequency [MHz] 



52 54 56 58 
Flux Coil Bias [jx A] 



FIG. 5: a) Points represent the observed linewidth of blue 
thermal sidebands as the JPA is tuned, employing a A' e = 
2tt x 13 kHz and n c = 1.6 x 10 6 coupling tone. The line is a 
theoretical prediction with no free parameters, b) From the 
same experiment, the observed changes in the oscillator center 
frequency and theoretical prediction are shown, c) Cooling 
efficiency (CE), as inferred from ring-down measurements and 
measured red and blue sideband powers. Left, CE for different 
coupling tone frequencies when JPA is far-detuned. Right, 
CE for A' 9 = 2tt x 13 kHz and different JPA tunings. Lines 
are theoretical predictions. Numbered call outs in right plot 
correspond to JPA tunings used in Fig. [6k. 



K n is the total network EM linewidth. In Fig. [5]3, left 
plot, we plot CE as a function of coupling tone fre- 
quency and with the JPA far-detuned. These are mea- 
sured by first inducing coherent oscillations in the me- 
chanics well above their thermal occupation through am- 
plitude modulation of the coupling tone at frequency £7, 
and then stopping the modulation and measuring the 
red and blue sidebands emitted by the network as the 
mechanical state re-equilibrates. Data is averaged over 
25 trials. The power emitted by both sidebands, P re d 
and Pbiue-, decays exponentially ("rings down") in time 
at rate T, while CE=(P Uue - P r ed)/{Pblue + Pred) is 
constant. The extremum of CE are ±0.42, implying a 
ft n /27r = 3 MHz-linewidth network coupled to the signif- 
icantly slower ^/2tt = 714 kHz oscillator, an unresolved- 
sideband network. Underlying this data is the expected 
CE versus coupling frequency, using independent network 
and oscillator calibrations. Driving the network as in 
Figs. [5^i&b and tuning the JPA, but now making ring 
down measurements, CE peaks at 0.98 at 56.4 /iA. A 
CE of 0.98 is only possible for networks with at most 
ft n /27r « 300 kHz, a resolved-sideband network and a 
10-fold reduction in n n consistent with the thermal data 
estimate above. Such a reduction was apparent from mi- 
crowave S\\ measurements in Figs. [4fsz[3]3, but CE repre- 
sents an effective Si 2 measurement from the mechanical 
bath port to microwave output (i.e. S\\ measurements 
test S MoM . [s] while CE, £12 measurements test S /ioT / l . [s]). 
CE also reaches a minimum of -0.04 at 57.8 /iA, indi- 
cating that the network resonance can also dip slightly 
below its far-detuned JPA resonance. An independently 
calibrated control systems prediction using 



CE = 



p{/i Th 



I 



2 i„ 



2 



(14) 



underlies the data. 



Finally, we demonstrate that the state of the JPA, and 
thus the control it effects, may be switched far faster 
than the oscillator can equilibrate. While above we 
controlled the JPA state using static flux biasing, be- 
low we apply both static flux biasing and an additional, 
A f e = 2tt x 10 MHz microwave tone, whose amplitude 
is dynamically switched and is not cancelled at the di- 
rectional coupler. At its strongest, this tone is too weak 
and off-resonant with the EMC to effect the oscillator, 
but is sufficiently strong and resonant with the wider- 
band TKC to strongly shift its center frequency. The 
rate with which the JPA may be switched in this manner 
should lie between the TKC linewidth (50 MHz) and the 
network's EM linewidth as a whole (>300 kHz ), orders of 
magnitude faster than the mechanical response rate (0.81 
Hz intrinsic linewidth). By performing similar ring down 
measurements starting from an even larger mechanical 
coherent amplitude (inducing ~100 kHz pk-pk resonance 
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shift in the EMC), this switched control may be observed 
with high visibility and bandwidth. In Fig.|6^, the power 
measured in the blue and red mechanical sidebands is de- 
picted in time where, from to 25 ms, the network is in 
a resolved-sideband and blue detuned state (state 1 in 
Fig. J6Jd, with CE called out in Fig. [5]^), while from 25 
to 40 ms, the network is in an unresolved-sideband and 
slightly red detuned state (state 2). In the first segment, 
almost all the mechanically-scattered power is in the blue 
sideband, and the total power emitted decays in time, 
indicating oscillator cooling (the large amplitude motion 
employed here induces a variation in the EMC center fre- 
quency of order the state 1 linewidth. As a consequence, 
the ring down power is more linear than exponential in 
time). In the second segment, the two powers are roughly 
equal, with the red sidband slightly greater, and the to- 
tal power hardly decays. The inset depicts the transition 
between states 1 & 2, through which the blue sideband 
stays constant, but the red sideband power jumps by an 
order of magnitude. This is only possible because state 
2 has a much larger linewidth than 1 (Fig. ^p). The 
transition occurs at least as fast as the 25 kHz detection 
bandwidth, much faster than the mechanical decay rate 
in state 1 (r <C T ~ 2tt x 15 Hz) or 2 (r > T = 2tt x 0.1 
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FIG. 6: Dynamic JPA-control. a) Red and blue sideband 
powers ringing down from a very high coherent amplitude. 
From 0-25 ms (25-40 ms) the network is in state 1 (state 
2), see Figs. |5j; & [6j). Inset, detail of network transition, 
depicted with 25 kHz bandwidth, b) Microwave probe Su 
measurements for the network states employed in a, with ar- 
rows indicating coupling tone and sidebands, c) Switched re- 
thermalization of a cooled oscillator. At 100 ms, JPA is tuned 
from T ^> To (cooling) state to F = To (bath-thermalizing) 
state. At 600 ms, JPA is switched to the cooling state for 
6 ms only. Plot represents the ensemble distribution of blue 
sideband amplitudes in time for 300 trials, d) Sideband am- 
plitude variance for 6 ms=2/r and 3 ms=l/r cooling- interval 
experiment ensembles. 



Hz). 

Switched control over incoherent thermal motion may 
also be demonstrated, at the expense of reducing the de- 
tection bandwidth to 100 Hz (so that 10 min averaging 
times yield reasonable signal to noise ratios). In Fig.[6]>d, 
the coupling tone and JPA are tuned such that the me- 
chanics are cooled for several seconds with T = 2tt x 56 Hz 
(the most aggressive, resolved-sideband cooling configu- 
ration in Fig. [5^). At time 100 ms, the JPA is rapidly 
far-detuned from the EMC, so that the network is on- 
resonance with the coupling tone and is in an unresolved 
sideband limit. For 500 ms, the mechanical state then 
thermalizes with its bath at rate IV After 500 ms, the 
JPA is briefly tuned to the cooling state for 6 ms, and 
then back to the thermalization state for a final 400 ms. 
Fig. [6)3 depicts the 300-trial ensemble distribution of con- 
tinuously monitored blue sideband amplitudes. Before 
100 ms, the distribution is dominated by amplifier noise, 
but broadens by a factor of 3 during the thermalization 
periods. The sideband amplitude is proportional to me- 
chanical position (by a factor that differs for the cooling 
and thermalization states) , and the broadening of the dis- 
tribution is indicative of a warming oscillator (to ~ 750 
phonon occupation, or 26 mK, or (63 Hz) 2 variance in 
the EMC resonance, although precise knowledge of the 
thermal occupation is incidental to our purposes). The 
position variance of the hot r = 2tt x 56 Hz oscillator is 
expected to drop by e -2 in 6 ms, by e _1 in 3 ms, and the 
variance of the cold, To-linewidth oscillator should rise to 
1 — e _1 of the bath equilibrium value in 200 ms. Fig. [6]i 
depicts the ensemble variance of 6 & 3 ms cooling-interval 
experiments and is consistent with these expectations. 



CONCLUSION 

We have demonstrated a small coherent feedback net- 
work of superconducting microwave devices that provides 
a degree of dynamically flexible control over a mechan- 
ical oscillator previously unavailable to the field of elec- 
tromechanics. The network's operation has two natu- 
ral interpretations: a dynamically-tunable input coupler 
[42] , or a gain-1 amplifier that receives and feeds coherent 
signals back to the oscillator and readout channel with 
a controlled phase shift [9]. Despite the simplicity of 
the components and construction (SMA cabling between 
pre-existing superconducting devices) , the network is too 
complex to be modeled using traditional electromechan- 
ical techniques (e.g., such approaches assume a single 
resonator mode, while our network features coupled res- 
onators) [33-35 . However, it is efficiently and intuitively 
modeled using a linear coherent feedback approach [5]- 
|9]. Although we only demonstrate this network's opera- 
tion in the classical regime of mechanical motion and EM 
fields, on the basis of the well-known connections between 
classical and quantum dynamics in linear coherent sys- 
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terns [5H21 E] , we expect that the essential mechanisms 
of the network should work analogously in the presence 
of unambiguously quantum fields and states. We have 
not yet probed the quantum regime simply because of 
our use of a HEMT (i.e. non-quantum limited) amplifier 
for readout and the dynamical richness to be explored 
first in the classical regime. 

There are several worthwhile avenues for future inves- 
tigations. First, the TKC may be pumped by a third mi- 
crowave tone to add a parametric gain-and-squeezing el- 
ement [25] (or a more general Kerr nonlinearity [13] ) into 
the feedback and read-out dynamics. Recent theoretical 
work suggests that coherent feedback to an electrome- 
chanical system from a parametric gain controller can 
outperform any type of ideal measurement-based feed- 
back or passive coherent controller [7J [9]. Second, the 
ability to dynamically and continuously modulate the 
network's coupling to its I/O port could be leveraged 
to shape the waveform of information-carrying signals 
read into and out of the EMC. This capability could fa- 
cilitate high-fidelity coherent state transfer between the 
oscillator and arbitrary coherent devices either "up-" or 
"down-stream" from the network [4TH44] . And finally, 
the model could be considered from the perspective of 
well-developed theories of classical optimal and robust 
control [TJ |6j [9] . These sorts of investigations are likely 
to yield recommendations for a more precisely-controlled 
network construction. In particular, different applica- 
tions are helped and hindered by different round trip 
phase shifts over the EMC and TKC network branches. 
While these phases were uncontrolled by the use of bulky 
SMA interconnections in our work above, a more inte- 
grated network could be constructed with much better 
precision. 
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